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Abstract 



: 

We consider scenarios with large flavour violating entries in the squark mass matrices focusing on the 
mixing between second and third generation squarks. These entries govern both, flavour violating low 
energy observables on the one hand and squark and gluino decays on the other hand. We first discuss 



I the constraints on the parameter space due to the recent data on B mesons from the B factories and 
' Tevatron. We then consider flavour violating squark and gluino decays and show that they can still 
■ be typically of order 10% despite the stringent constraints from low energy data. Finally we briefly 
^ . comment on the impact for searches and parameter determinations at future collider experiments 
O I such as the upcoming LHC or a future International Linear Collider. 
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1 Introduction 



Theoretical arguments like the hierarchy problem motivate the general expectation that the ex- 
periments at the LHC will lead to discoveries of new degrees of freedom at the TeV energy scale. 
The precise nature of this new physics is unknown, but it most probably will answer some of the 
fundamental questions related to the origin of electroweak symmetry breaking. 

Rare B and kaon decays (for a review see [1,2]) representing loop-induced processes are highly 
sensitive probes for new degrees of freedom beyond the SM establishing an alternative way to search 
for new physics. However, this indirect search for new physics signatures within flavour physics takes 
place today in complete darkness, given that we presently have no direct evidence of new particles 
beyond the Standard Model (SM). But the day the existence of new degrees of freedom is established 
by the Large Hadron Collider (LHC), the searches for anomalous phenomena in the flavour sector 
will become mandatory. The problem then will no longer be to discover new physics, but to measure 
its (flavour) properties. 

Thus, within the next decade an important interplay of flavour and high-p^ physics most probably 
will take place. For example, within supersymmetric extensions of the SM, the measurement of the 
flavour structure is directly linked to the crucial question of the supersymmetry-breaking mechanism 
as the soft SUSY breaking terms are the source of flavour structures beyond the SM. LHC has 
the potential to discover strongly interacting supersymmetric particles up to a scale of 2 TeV and to 
measure several of their properties [3-6] . This information can be used for a refined analysis of flavour 
physics observables indicating possible flavour structures and, thus, give important information for 
distinguishing between models of supersymmetry breaking. 

Data from K and Bd physics show that new sources of flavour violation in s ^ d and h ^ d are 
strongly constrained, while the possibility of sizable new contributions io h ^ s remains open. We 
also have hints from model building: flavour models are not very effective in constraining the 6 — > s 
sector [7]. Moreover, in SUSY-GUTs the large mixing angle in the neutrino sector relates to large 
mixing in the right-handed h-s sector [8-10]. 

As we explicitly show in this paper, such flavour information on the h ^ s observables is comple- 
mentary to the high-pT data of the LHC. In fact, squark and gluino decays are governed by the same 
mixing matrices as the contributions to flavour violating loop transitions. This allows for possible 
direct correlations between flavour non-diagonal observables in B and high-p^- physics. We already 
anticipate that the present bounds on squark mixing, induced by the low-energy data on 6 — ^ s 
transitions, still allow for large contributions to flavour violating squark decays at tree level. Due 
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to the restrictions in flavour tagging at the LHC, additional information from future flavour experi- 
ments will be necessary to interpret those LHC data properly. Also the measurement of correlations 
between various squark decay modes at a future ILC would provide information about the flavour 
violating parameters. 

Flavour violating squark and gluino decays have already been considered some time ago [11]. In 
the present paper we generalize and update the previous analysis including additional experimental 
constraints. More recently, similar work on the charged Higgs boson production was presented; it 
was shown that squark mixing can signiflcantly change the production pattern [12]. In another work 
the interesting question was addressed if high-p^- data can contribute to the solution of the flavour 
problem [13, 14]. Recently, it was shown that one can also derive signiflcant bounds on flavour- 
violating parameters of the squark sector by requiring that the radiative corrections to the CKM 
elements do not exceed the experimental values [15]. Flavour violating squark production has been 
studied in [16-18] where it has been found that flavour violating production can be sizable after 
taking into account constraints from the b ^ s observables and the anomalous magnetic moment of 
the muon. 

The paper is organized as follows: In Chapter 2 we review the necessary information about flavour 
mixing in the MSSM and about decays of squarks and gluinos in order to introduce our notation. 
In Chapter 3 we present our phenomenological results and in Chapter 4 we discuss their impact for 
LHC and ILC. 

2 Preliminaries 

2.1 Decays of squarks and gluinos 

In the study of squark decays two scenarios can be distinguished depending on the SUSY spectrum: 
1. rrig > TTiq. [q = d,U] i = 1, . . . ,Q): In this case the gluino will mainly decay according to 

g^djdi, g-^UjUi (1) 
with dj = {d, s, b) and Uj = {u, c, t) followed by squark decays into neutralino and charginos 

Ujxl , djxi , di djxl , UjxY ■ (2) 
In addition there can be decays into gauge and Higgs bosons if kinematically allowed: 

Ui^Zuk, H^Uk 1 W~^dj , H^dj] di^Zdk, H^dk, W~Uj , H~Uj (3) 
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where = ^4°), k < i, j = 1, . . . ,6. Note, that due to the fact that there is left- 

right mixing in the sfermion mixing, one has flavour changing neutral decays into Z-bosons at 
tree- level. 



2. nig < niq. (g = d,u; i 



1, . . . , 6): In this case the squarks decay mainly into a gluino. 



(4) 



and the gluino decays via three-body decays and loop-induced two-body decays into charginos 
and neutralinos 



with i,j = 1, 2, 3, / = 1, 2 and = 1, 2, 3, 4. Note that the first two decay modes may contain 
states with quarks of different generations of quarks. 

Explicit formulas for the partial widths including flavour effects can be found in Ref. [17]. The flavour 
mixing final states of the decays listed above are constrained by the fact that all observed phenomena 
in rare meson decays are consistent with the SM predictions. As we will show in the next sections, 
there are regions in the parameter space where the flavour violating decay modes can be even of the 
order of 10%. 

2.2 Flavour changing neutral currents 

Within the Minimal Supersymmetric Standard Model (MSSM) there are two new sources of flavour 
changing neutral currents (FCNC), namely new contributions which are induced through the quark 
mixing as in the SM and generic supersymmetric contributions through the squark mixing. The 
latter is described by their mass matrices. 



where / stands for up- or down-type squarks. In the super-CKM basis in which the quark mass 
matrix is diagonal they read as 



g dj di xl , Uj Uixl, Q^u^dixf , 9^9 xl 



(5) 





DfLL = (Tsj - e/ sin^ 9w) cos(2/3)m| 



z ■> 



DfRR = e/ sin^ 9w cos(2/3)m| 



(7) 
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for the D-terms, 

FfLL,ij = FfRR^ij = mjSij , FfRL,ij = - firuiSij (tan py'^'^-'-f (8) 

where rrii are the corresponding quark masses, Cf the electric charge, and T-^j the weak isospin of 
the corresponding left-handed quark. In this basis the F- and D-terms are flavour diagonal and all 
flavour violation beyond the CKM resides in the soft SUSY breaking terms: 

MIll = V^Kj^Ml^^VcKM = ml^VjmlV,, (9) 
MlnR = m\^V\mf\J,, M^j^j, = ml ^ Ulmf , (10) 

MIlr = vjV2fD^v,/V2UlTEV,, MlL^ = V2/V2fu^V2/V2UlT^Vu, (11) 

where the un-hatted mass matrices mg ^ ^ and trilinear interaction matrices T^.zj are given in the 
electroweak basis. The transformations Vu^d and Uu^d just bring the quarks from the interaction 
eigenstate basis to their mass eigenstate basis, so Vckm = Vu'^d- The relation between Mu,ll and 
Md^LL is due to SU (2) l gauge invariance. The T-matrices are in general non-hermitian|*] 

These additional flavour structures induce flavour violating couplings to the charginos, neutralinos 
and gluinos in the mass eigenbasis, which give rise to additional contributions to observables in the 
K and B meson sector. The low-energy observables can thus be used to constrain the size of the 
off-diagonal elements of the mass matrices Mj and Mj i^^. Some of the flavour- violating 

terms may turn out to be poorly constrained. Thus, it is suitable to rely on the mass eigenstate 
formalism, which remains valid - in contrast to the mass insertion approximation - even when the 
intergenerational mixing elements are large. The diagonalization of the two 6x6 squark mass 
matrices squared and M."^ yields the eigenvalues and [k = 1, 6). The corresponding 
mixing matrices R{j relate the mass eigenstates with the electroweak eigenstates 

Uk = Rl.uY , 4 = R%dY (12) 

where uf G {ml, cl, ^l, ur, cr, Ir} and d'j'" G {di, sl, h, dR, sr, Ir}, respectively. 

As usual, the flavour off-diagonal elements of the squark mass matrices are normalized by the 
diagonal ones. One uses the average of the diagonal elements (trace of the mass matrix divided by 
six) in the up and down sector, denoted by m?. The observables can then be studied as a function 
of the normalized off-diagonal elements 

dLL,ii = -2 , = . K^T3) (13) 



*We follow here the conventions of the SUSY Las Houches Accord, for more details see Ref. [19]. 
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where / is either u oi d for w-squarks and (i-squarks, respectively. We emphasize that a consistent 
analysis of the bounds should also include interference effects between the various contributions, 
namely the interplay between the various sources of flavour violation and the interference effects of 
SM and various new-physics contributions, (see Ref. [20]). 

At present, new physics contributions to s ^ d and h ^ d transitions are strongly constrained. In 
particular, the transitions between first- and second-generation quarks, namely FCNC processes in 
the K system, lead to very strong constraints on the parameter space of various new physics models. 
However, most of the phenomena involving 6 — > s transitions are still largely unexplored and leave 
open the possibility of large new physics effects, in spite of the strong bounds of the three most 
important b s observables, namely, the inclusive decay modes B Xg'y and B Xsi~^i~ and the 
Bs - Bg mixing. It is well-known that the gluino contributions to these two inclusive decay modes 
are mostly sensitive to 5d,2z,LR and 5d,2z,RL in the leading order approximation [21-23]. However, 
the restriction to a single 5 parameter and to the gluino contribution lead to unrealistically strong 
constraints as was already pointed out in Ref. [20]. There are also sensitivity to combinations of 
other 5 parameters and, moreover, depending on the precise choice of the supersymmetric flavour- 
diagonal parameters, other supersymmetric contributions which are also sensitive up-squark mixing 
parameters. The gluino contributions to Bs — Bg mixing are mainly sensitive to the combination 
^LL^d,RR and 6d,LR^d,RL- The latter is not really relevant due to the constraints already induced 
by the inclusive decay modes. But also here sensitivity of that observable to other squark mixing 
parameters are not negligible. 

Moreover, there is an impact of squark mixing including the third generation on the lightest Higgs 
mass and the p parameter which also includes the 5 parameter in the up-squark sector (see for 
example [24]). There are the Tevatron bounds on squark masses and certain constraints from dark 
matter phenomenology. Also the diagonal T entries (see Eq. [TTl) have to be tuned that they lead to 
perturbatively quark masses. Finally, we recall that there are also constraints on some 5 parameters 
from theoretical considerations, namely that one should avoid colour and charge breaking minima [25]. 
All resulting numerical bounds and constraints will be listed below. 
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3 Phenomenological analysis 



3.1 Benchmark points, experimental and theoretical constraints 

In our present analysis, we first fix the flavour-diagonal set of parameters and then we vary the 
flavour-nondiagonal parameters and explore the bounds on those parameters by theoretical and 
experimental constraints. We fix the flavour-diagonal parameter set following three very popular 
SUSY benchmark points, namely SPSla' [26] ,/" and 7 [27]. The first two points are mSUGRA 
points whereas the third one is a mSUGRA point with non-universal Higgs-mass parameters at the 
GUT-scale. SPSla' contains the lightest spectrum with squarks around 500 GeV and rrig around 
600 GeV, tan/3 = 10, followed by 7 with squark masses around 600 GeV, rrig around 580 GeV, 
tan/5 = 20, and /" with squark masses around 730 GeV, rrig around 1000 GeV, tan/5 = 35. For 
more information on the benchmark points see Refs. [26,27]. All of them are consistent with WMAP 
data [28] and measurements of the anomalous magnetic moment of the muon. 

On the flavour-nondiagonal parameter set we pose the following concrete constraints: 

• We explicitly check if all our data points fulfill the theoretical vacuum stability constraints [25]. 
Actually, all our data points fulfill the conservative condition that in the LR and RL submatrix 
off-diagonal elements are not larger than diagonal elements, which implies the usual vacuum 
stability conditions [25]. 

• Among the constraints from electroweak precision data, the most important one is rriho > 114.4 
GeV, where we add 3 GeV to the theoretical prediction of rriho as a measure of the theoretical 
uncertainty [29-31]. Furthermore we require that = 80.40 ± 0.03 GeV where we take as 
input mz = 91.187 GeV, Gp = 1.16639 ■ 10"^ GeV^ and aenXO) = 1/137.0359895 [32]. 

• Squark Tevatron bounds on squarks are of the order of 250 GeV depending on the SUSY 
spectrum [33]. 

• The explicit experimental constraints from the most important flavour observables we use in 
our analysis are [34-39]: 



2.67 < 



13.5 < 



Br{B Xs-f) X 10^ 
AMbs ps 



< 21.1 



< 4.29 



(15) 
(16) 



1.05 < 



,2 X 10' 



< 2.15 



(17) 



BR{Bs /i+/i~) X 10' 



|8 



< 5.8 



(18) 



7 



Those bounds include experimental and theoretical errors which are linearly added. Explicitly 
our bounds are the experimental 95% bounds where twice the SM error is added in order to 
take into account uncertainties of the new physics contributions in a conservative way. We have 
also checked that the recent experimental data on B ^ ru do not give additional constraints. 

• Another requirement is that the lightest stable particle (LSP) should be neutral (but not 
necessarily be the neutralino) which however is fulfilled once the experimental bounds on the 
squark masses are taken into account. 

For the numerical evaluation we use an updated version of SPheno [40] which has been extended 
to accept fiavour mixing entries in the sfermion mass matrices. The masses and mixings of all SUSY 
particles are calculated at the one-loop level where the formulas of ref. [41] have been extended for the 
fiavour- mixing case [42]. The masses of the neutral Higgs bosons are calculated at the two- loop level. 
For the 5-physics observables we use the formulas of Refs. [34,43-45] for B ^^7, Refs. [39,44-46] 
for B — s> Xgl'^l" , and Ref. [38] for AM^^ and B^ /i'^/.i". The branching ratios for squark and 
gluino decays are calculated using tree-level accuracy but with running couplings evaluated at the 
scale Q = m where m is the mass of the decaying particle. 

3.2 Results 

3.2.1 Constraints on fiavour parameters 

In our analysis we want to identify the fiavour violating decay channels of squarks and the gluino with 
potentially large branching ratios. And we are interested how the various decays can be combined 
to obtain information on the parameters. 

In Table [T] we present the allowed ranges of the fiavour- diagonal parameters 6ij in the b-s sector 
using the so-called 'one- mass-insert ion approximation' where just one 6 parameter is assumed to 
be active and all the others are set to zero. Note that we have calculated observables using the 
exact diagonalization of the mass matrix of eq. Our findings are consistent with the results 
of previous analyses in Refs. [12,20,21,47,48], in particular the pattern of the constraints on the 
up-quark sector. Here it turns out that the parameter 6u,lr,23 receives bounds from the fiavour and 
electroweak constraints while the other two parameters 6u,lr.32 and 6u,rr,23 are unconstrained in 
scenarios with low and moderate tan /5 and weakly constrainted once tan j3 gets large. 

In Figure □ we present regions in the Sd,RR,23-SLL,23 (left) and in the 5d,LR,23-^d,LR,32 planes (right) 
consistent with experimental data for the SPSla' benchmark point. We show here the lines corre- 
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Table 1: Allowed ranges of the 6 parameters in the neighborhood of the benchmark points SPSla', 
7 and I" taking into account the experimental information given in the text and assuming that only 
one flavour-mixing parameter is present. The regions correspond to 95% CL. 





SPSla' 


7 


I" 




(-0.05,0.03) 


(-0.037,0.005) 


(-0.06,0.001) 


Sd,RR,23 


(-0.43,0.66) 


(-0.29,0.48) 


(-0.5,0.45) 


Su,RR,23 


(-0.7,0.7) 


(-0.54,0.43) 


(-0.55,0.45) 


Su,LR,23 


(-0.16,0.08) 


(-0.16,0.06) 


(-0.35,0.05) 


5u,LR,32 


(-0.7,0.54) 


(-0.5,0.2) 


(-0.7,0.27) 


Sd,LR,23 


(-0.0047,0.0046) 


(-0.006,0.001) 


(-0.01,0.0015) 


Sd,LR,32 


(-0.019,0.02) 


(-0.015,0.015) 


(-0.004,0.003) 



sponding to the constraints from b ^ s'-f using red (full) lines and AMb^ using magenta (dashed) 
ones. The areas consistent with all constraints correspond to the blue (dark) areas. In both cases the 
complete plane is compatible also with the constraints due to the observable b sl'^l" , Bg fi'^fi" 
and Bu — > tu. 

One sees the regions allowed by all constraints corresponding to the blue (dark) area are mainly 
along the axes as one would naively expect. Note, however, that the constraint on a parameter gets 
more involved once one allows for additional sources of flavour violation. For example. Figure [T^ 
shows that for small negative Sll,23 the constraint on 5d,RR,23 becomes first weaker; for larger (5ll.23 
the bound is stronger again. 

In both planes of Figure 1 we find an additional region consistent with the data in the upper right 
corner where both flavour mixing parameters are sizable. However, in these regions clearly some 
cancellations between various contributions to the observable take place. In particular the gluino 
and the chargino contributions for AMb^ are large and of opposite sign. Their modulus can reach 
about twice the size of the SM contribution. 

The other two benchmark points show very similar features but for the fact that in case of /" also 
Bs fi^fi" would exclude part of the allowed regions. 
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a) b) 




-0,10 -O.OS OOO O.OS 0.10 -0,015 -0,010 -OOOS 0,000 0,005 0010 OOlS 



^LL,23 ^D,LR,23 

Figure 1: Allowed regions a) in the Sd,RR,23-SLL,23 and b) in the Sd,lr,23-Sd,lr,32 plane for SPSla'. 
All other than the shown flavour off-diagonal elements are put to zero. The lines correspond to: full 
red lines 6 — S7 , 2.6710"^ and , b ^ s'j , 4.2910""^, dashed magenta lines lAjv/^J = 13.5 ps^\ and 
|Ajv/bJ = 21.1 ps^^. The blue (dark) area shows the regions consistent with all data at 95% CL. 

3.2.2 Flavour violating decays of squarks and gluinos 

In the following we discuss the effect of the flavour mixing parameters on the decay properties of 
squarks and gluinos. For the discussion of the basic features we will first take two points in the 
allowed region of Figure [T^). Afterwards we discuss various parameter dependencies. Note that we 
often show only part of the kinematically accessible final states and, thus, the branching ratios shown 
do not sum necessarily up to 1. For the flavour diagonal entries we will use the point SPSla' but 
again similar features are found for the other benchmark points as well. The main difference is due 
to the kinematics. 

The two study points chosen are characterized by Sll,23 = 0.01 and Sd,rr23 = 0.1 (point I) and 
^LL,23 = 0.04 and 5d,rr23 = 0.45 (point II) respectively. Study point II is characterized by large 
cancellations of the SUSY contributions to 5-physics observables. In Table H] we give a summary of 
the various branching ratios and in Table [3] we display the masses. For comparison also the masses 
without flavour mixing parameters are given. 

The flrst feature to note is that the bounds on 5ll are already so strong that there are small 



10 



Table 2: Branching ratios larger than 1% for two study points. The flavour diagonal entries are 
according to SPSla'. Ui decays are like in SPSla' [26] and in both scenarios BR(c/3 — > xld) = 99.1%. 



decaying 




final states and corresponding branching ratios in % for. 




J- • 1 

particle 


I. 5 


LL, 23 = 


= 0m,6D,RB23 = 0.1 


If- ^LL,23 = 


= 0.04, 6d,rr23 = 0.45 


di 


X\h, 4.4 


~ n I. on 

X20, 29.8 


Xit, 37.0 


36.8 


Xib, 42.2 


~0i inn 

X20, 10.9 




uiW , 


27.7 






Xit, 9.6 






3 
"2 




8.0 


~0l. C A 

Xib, 6.4 


~0 inn 

X20, 19.0 


2.1 


X20, 27.3 


Xit, 34.6 






1.1 


Xlb , 1.8 


xrt, 24.6 


uiW , 33.2 








uiW- , 


38.9 












d^ -H> 


XiS 


9.1 


Xib, 6.3 


X2S, 25.3 


2.3 


Xarf, 31.7 


Xi-u, 59.7 






2.1 


Xic, 47.3 


" T T 7" — AO 

uiW , 4.8 


~ — on 

Xic, 3.0 


X2U, 2.3 




«5 


Xid 


2.3 


X^d, 31.7 


Xiu, 59.9 


Xis, 2.2 


X^s, 30.7 


Xiu, 2.9 




XiC 


2.8 


X2U, 2.3 




Xfc, 58.5 


X2C, 2.3 




da 


X\s 


3.1 


X^s, 30.6 


xr^, 2.7 


19.7 


Xib, 18.8 


2.9 




XiC, 


58.1 


X2C, 2.4 




xlb, 2.9 
^6, 39.8 


X2-t, 5.8 
uiW~ , 5.5 


^s, 2.2 


9 


Uit, 


19.2 


U2C, 8.2 


u^u , 8.3 


■uit , 13.5 


'U2C, 5.8 


usu , 5.8 






4.2 


U5C, 4.2 




U4C, 2.6 


MsM, 2.6 






dis 


1.4 


dib, 20.6 




Jis, 21.1 


dib, 22.7 






d2S 


6.3 


d2b, 9.0 


(isfi, 8.3 


d2b, 14.0 




c/sc?, 5.9 




d^s 


2.3 


d^b , 1.3 


Jes, 2.8 


(i4(i, 2.3 


(isrf, 3.3 





effects only on masses or branching ratios due to this flavour-mixing parameterjH This implies that 
in our examples the masses and branching ratios of the u-type squarks are hardly altered compared 
to the SPSla' point. Things are different for the ci-type squarks and consequently also for the gluino 
as the 5d,_Rij,23 parameter can be sizable. 

The relative size of the branching ratios in Table [2] can be understood by the nature of the various 
squarks mass eigenstates. In point I one finds di ^ 6/, with a small admixture of b^, d2 — b^ with 
small admixtures of sr and bi, d^ ^ dR, d^ ^ sr with admixtures of sl and b^, d^ c:=l di and d^ c:=^ sl 
with a small admixture of sr. Thus, larger flavour effects are visible in the decays of ^2 and ^4 where 

t Note that this strong constraint would be weakened if more than two deltas are active as shown in the study 
point given in ref. [11] that contains such an example. 
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Figure 2: Composition of a) di=i and b) di=2 as a function of 5d,RR,23i the flavour diagonal parameters 
are the ones of SPSla'. The flavour basis is given by (c/l, Sl, hi, (Ir, sr, The hues correspond to 
the final states: dashed green line to iW^r P, dotted red line to l-Rf ; P, and full black line to iR'^r P. 



the flavour violating decay branchings ratios ^2 Xi^ and Xi^ are of the order of 10%. This 

structure is also the reason for the relative importance of the flavour violating decays of the gluino. 
As a side remark we note that the flavour violating decays of the first generation squarks and of the 
2nd/3rd generation squarks into the first generation quarks are due to CKM quark mixing. 

In point II the situation is more complicated due to the larger flavour mixing parameters. With 
respect to the nature of the d-tjpe squarks we find that di and rfg are strongly mixed states consisting 
mainly of sr and Br with a small admixture of whereas the other states are mainly electroweak 
eigenstates: d2 — b^, ds ^ dR, d^ ~ di and d^ ^ sl. In this scenario the flavour violating final states 
can even reach about 40% in case of d2 — > x^s and about 20% for rfg — > Xi^- The differences for the 
gluino decays between these two points is not only due to the different mixing in the (i-squark sector 
but also due to the different kinematics as can be seen from Table El 

The effect of the flavour mixing parameters on gluino decays is strong except for 5d,LR,ij- We 



Table 3: Squark masses in GeV for SPSla' and our two points under study. The flavour diagonal 
entries are according to SPSla'. Note that rriu^ ~ m^^g and rriu^ ~ mj^g. 















mi 




mi 
de 












SPSla' 




506 


546 


547 


547 


570 


570 


367 


547 


565 


586 




= 0.01,5d,/jij,23 


= 0.1 


503 


525 


547 


569 


570 


570 


366 


547 


565 


586 


n. 5lL,23 


= 0.04, 5d,RR,23 


= 0.45 


422 


509 


547 


570 


572 


641 


366 


547 


565 


587 
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Figure 3: g decays as a function of 5(1,rr,23i in a) and b) all other 5s are zero whereas in c) and d) 
^u,RR,23 = 0.2 in addition but all other 5 parameter are zero. The lines in a) and c) correspond to 
the final states uic (long dashed blue line), U2C (short dashed green line), dis (dotted red line), and 
(full black line); the lines in b) and d) correspond to the final states uit (long dashed blue line), 
(short dashed green line), dih (dotted red line), and (full black line). 



first investigate how the nature of the various squark states change when varying the flavour mixing 
entries. In Figure [2] we display the composition of the two lightest (i-squarks as a function of 5d,_RR,23- 
The lightest state is mainly a hi state for \5d,RR,23\ ^ 0.15. At \5d,RR,23\ ^ 0.2 it is a strongly mixed 
state consisting of 6^, 6r and sr whereas for larger values of \5d,RR,2-i\ it is mainly an admixture of 
hfi and sr. The second state is mainly a h^i state for |5d,_Ri?,23| ^ 0.05. The kinks around 0.05 stem 
from the fact that the flavour mixing is about as large as the left-right mixing in the sbottom sector. 
Around \5d,RR,23\ ^ 0.15 the hi component starts to dominate; for \5d,RR,23\ ^ 0.3 this state is almost 
a pure hi. 

This behaviour is reflected in the gluino decays as can be seen in fig. O where we show the gluino 
decay branching ratios into the two lightest (i-squarks and into the two lightest w-squarks as a 
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Figure 4: Composition of a) Ui=\ and b) Ui=2 as a function of 5u,_Ri?,23, the flavour diagonal parameters 
are the ones of SPSla'. The flavour basis is given by {uL,CL,iL,UR,cji,iji). The lines correspond to 
the following final states: dashed green line to |i?"r P, dotted red line to \Rf 



, . 1^, and full black line 



to p. 



function of 5d,R_R,23- As already indicated in the discussion above, the dependence is twofold. The 
larger the off-diagonal elements are, the larger is the corresponding flavour violating coupling; but 
also kinematical effects play an important role, e.g. the larger 5d,RR,23, is, the lighter di gets leading to 
an increased phase space. The kinks and the crossing of the various lines correspond exactly to the 
changes in the nature of the d-squarks. Figure [St shows in addition that there are scenarios where 
one can have sizeable flavour violating decays into rf-squarks and into w-squarks at the same time. 

Similar features can be observed for the case of large flavour mixing in the u-squark sector. In 
Figure m we display the flavour content of the two lightest w-squarks as a function of 5u,rr,23,- There 
are qualitative differences compared to the (i-squarks because (i) the left-right mixing in the top- 
squark sector is signiflcantly large and (ii) the lightest state Ui in the flavour conserving case is a 
mixture of t/j and and U2 — cr- We see that with increasing (5„,_rr,23 not only the mixing between 
Ir and cr get larger as expected but that at the same time the admixture of II decreases quickly. 
The second state is not, as one would naively expect, a strongly mixed state consisting of Ir and cr 
but a mixture between and cr and only a negligible admixture oUr. The reason for this surprising 
feature is that the left-right mixing in the top-squark sector is larger than 5n,R_R,23 for the complete 
range shown. 

Again this is reflected in the gluino decays shown in flg. [3] where we display the gluino decay 
branching ratios as a function of 5u,r_r,23- Here the main flavour effect is the increase of the branching 
ratio into uic (long dashed blue line on the left side). The expected increase of U2t does not take 
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Figure 5: g decays as a function of Su,rr,23 and all other 6 parameters are zero. The lines in a) 
correspond to the final states uic (long dashed blue line), U2C (short dashed green), and (full 
black); the lines in b) correspond to the final states Uit (long dashed), (short dashed green), dib 
(dotted red), and (full black) 



place due to kinematics as + > rrig for most of the range shown. One would expect that 
the branching ratio for Uit should decrease for increasing |5u,_Ri?,23| as the top-squark components 
decrease. However, this decrease is overcompensated by the increased kinematics as the mass of ui 
drops to a value about 200 GeV. 

The squarks stemming from the gluino decays will decay further and in Figures Eland [7] we display 
the corresponding branching ratios of the lightest M-squark and (i-squark, respectively. The behaviour 
of ui decays in Figure [6] is a direct consequence of the above described change of its nature with 
increasing \6u,rr,23\- But also its mass decrease again implies that if \Su,rr,23\ increases first the final 
state and then also the final states xfb and Xit g^t kinematically forbidden. The kinks in the 
various lines are a direct consequence of these kinematical effects. Finally, for |5t(,i?_R,23| ^ 0.6 only 
the final state XiC is possible. We have checked that three and four body decay modes, such as 
bWxi [49], lub [50,51] or bff'xi [52], are suppressed and have branching ratios below 1%. 

In Figure [7] we display the branching ratios of di. The only sizeable flavour violating decay channel 
is Xi^ whereas XiC is suppressed because we consider here only flavour mixing between the -R-squarks. 
The above mentioned level crossing at |5d,ij_R,23| — 0.15 is the main reason for the drop of the final 
states X2^ and Xi^ besides kinematical effects. 
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Figure 6: Ui decays as a function of 6u,rr,23, the flavour diagonal parameters are the ones of SPSla'. 
The hues in a) correspond to the flnal states (dashed green hne), (dotted red), and xfs (full 
black); the lines in b) correspond to the flnal states x^t (dashed green), (dotted red), and xfb. 
(full black). 



4 Impact for LHC and ILC 

We have seen in Figure 1 that constraints due to low-energy observables lead to interesting structures 
in the parameter space and the only certain parameter combinations are allowed in the case one varies 
two parameters independently. Of course things are much more complicated once all off-diagonal 
elements are allowed to be non-zero at the same time. In the end the crucial question is how 
information on flavour violating squark and gluino decays can be used to pin down the underlying 
parameters. The answer to this question requires detailed information on masses, production cross 
sections and branching ratios. 

We have seen that squarks and gluinos can have flavour changing decay modes of typically O(10)% 
branching ratio. This clearly has an impact on the discovery strategy of squarks and gluinos as well 
as on the measurement of the underlying parameters at the LHC and a future international linear 
collider (ILC). For example, in mSUGRA points without flavour mixing one flnds usually that the 
left-squarks of the flrst two generations as well as the right squarks have similar masses. Large flavour 
mixing implies that there is a considerable mass splitting as can be seen in Table El Therefore, the 
assumption of almost degenerate masses should be reconsidered if sizable flavour changing decays 
are discovered in squark and gluino decays. 

An important part of the decay chains considered for SPSla' and nearby points are g bbj bbx^ 
which are used to determine the gluino mass as well as the sbottom masses or at least their average 
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Figure 7: di decays as a function of of 6d.RR,23, the flavour diagonal parameters are the ones of 
SPSla'. The lines in a) correspond to the final states Xi^ (dashed green line), x^-^ (dotted red), and 
XiC (full black); the lines in b) correspond to the final states Xi^ (dashed green), X2^ (dotted red), 
and Xit (full black). 



value if these masses are close [53]. 

In the latter analysis the existence of two 6-jets has been assumed stemming from this decay 
chain. In this case the two contributing sbottoms would lead to two edges in the partial distribution 
d{BIi{g bbxi)/dmbb where rribb is the invariant mass of the two bottom quarks. As can be seen 
from Figure [8] there are scenarios where more squarks can contribute and consequently one finds 
a richer structure, e.g. three edges in the example shown corresponding to study point I. Such a 
structure is either a clear sign of flavour violation or the fact that the particle content of the MSSM 
needs to be extended. Moreover, also the differential distribution of the final state bsxi shows a 
similar structure where the edges occur at the same places as in the bb spectrum but with different 
relative heights. This gives a non-trivial cross-check on the hypothesis of sizeable flavour mixing. 
Clearly a detailed Monte Carlo study will be necessary to see with which precision one can extract 
information on these edges. However, such a study is well beyond the scope of this paper. Obvious 
difficulties will be combinatorics because in general two gluinos or a gluino together with a squark 
will be produced and, thus, there will be several jets stemming from light quarks. However, one could 
take final states where one gluino decays into d-tjpe squarks and the second into stops or c-squarks. 
In the second case effective charm tagging would be a crucial. 

Similar conclusions hold for the decay chains g —>■ cui —>■ ctxi ^^'^ 9 ^ tui ^ ctx? analyzed in [54] 
and also for the the variable M]^ defined in [55] which sums up final states containing tbxf. 

At an ILC the situation should be considerable easier: First one can tune the center of mass 
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a) 10^ d{BR{g ^ bbx^/drubb b) 10^ d{BR{g hsxTl/dm^, 
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rribb nibs 

Figure 8: DifFerential distributions (i(BR(^ hbx^/dmbb and d{SK{g bbxi)/dmbb as a function 
of rribb = \J {pb + Pb)^ (nibs) for point 1 defined in section I3.2.2[ In b) tlie sum over tlie charges is 
sliown: BR(^ bsxi) + BR(^ — > bsxi)- 

energy so that one studies in principle the states one after each other. Secondly, one can polarize 
both, electrons and positrons, [56, 57] and thus one influences the production rates for signal and 
background. However, the ILC will be limited to energies in the one TeV range and, thus, most likely 
only part of the squarks can be explored at this machine. However, already this partial information 
combined with LHC data can give useful information [58]. Additional information can then be 
obtained at a multi-TeV e~^e~ collider such as CLIC [59] . 

5 Conclusions 

Flavour-violating low- and high-energy observables are governed by the same parameters in super- 
symmetric models. A particular important question is whether the soft SUSY breaking parameters 
can have additional flavour structures beside the well-known CKM structure. In this paper, we have 
analysed flavour-violating squark and gluino decays in view of the present flavour data and have 
shown that they can be typically of order of 10% in the regions of parameter space where no or only 
moderate cancellations between different contributions to the low energy observables occur. If we 
allow for larger new physics contributions, e.g. the same order as the SM contributions, in the flavour 
observables, then even flavour-violating branching ratios of up to 40% are consistent with the present 
data. We have checked that this is a common feature for a couple of SUSY benchmark points like 
SPSla', 7, and /". We have explicitly derived the pattern of flavour-violating decay modes from the 
specific structure of the flavour-violating parameters including kinematical constraints. Finally, we 
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have briefly analyzed the direct consequences for the search of supersymmetric particles in specific 
examples. The full exploitation of the necessary modifications in the particle search calls for detailed 
Monte-Carlo analyses. 
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